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Overview 


On 17 March 2013 at 03:50:54 UTC, NASA detected a bright impact 
flash on the Moon caused by a meteoroid impacting the lunar surface. 

There was enhanced meteor activity in Earth’s atmosphere the same 
night from the Virginid Meteor Complex. 

The impact crater associated with the impact flash was found and 
imaged by Lunar Reconnaissance Orbiter (LRO). 

A known crater location provides “ground truth” for testing the 
geolocation technique. 

Luminous efficiency estimates can be made by combining flash and 
crater measurements. A sanity check of photometric procedures and 
crater scaling relations is also possible. 
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3 Nov 2008 
00:11:06.144 
100 ms 
m R = 7.7 
0.1 kg 

S. Taurid (27 km/s) 



22 Apr 2007 
03:12:24.372 
133 ms 
m R = 6.7 
0.08 kg 

Lyrid (49 km/s) 



I 
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Jacobs Atypical flash on 17 March 2013 £EB 


17 Mar 2013 
03:50:54.312 
1.03 s 
m R = 3.0 
16 kg 
Virginid 


Flash info 

Detected with two 
0.35 m telescopes 

Watec 209H2 Ult 
monochrome CCD 



cameras 

— Manual gain control 
— No integration 
- r = 0.45 

Interlaced 30 fps video 

Saturated — ► needed 
saturation correction! 
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Photometry 
performed using 
comparison stars 

(see Suggs et al. 2014) 


Peak R magnitude 
saturation correction 





Saturated 
Peak m R = 4.9 
UNDERESTIMATED! 

* 


CORRECTION: 


2D elliptical Gaussian fit 
to the unsaturated wings 

Peak m R = 3.0 ± 0.4 

^ + 39 

Luminous energy = 7.1 24 x 10° J 

(Similar results for 2D elliptical Moffat fit) 
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Increased lunar activity 
on 17 Mar 2013 
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Observed by A. Kingery & R.M. Suggs; detected by R.J. Suggs Obsen/ecnDy A. Kingery & R.M. Suggs; detected by R.J. Suggs 




Jacobs Meteor activity on 17 Mar 2013 




SOMN 


NASA network 


Googlfc eSfth 


19 fireballs were observed by NASA & SOMN 
all-sky meteor cameras on 17 Mar 2013 
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19 fireballs seen on 17 Mar 2013 



Geocentric meteor radiants color-coded by 
speed with a tight cluster of 5 with: 


Virginid Complex 
at X=356.6 



meteors 

NVI 1 

EVI 2 

«g(°) 

184.1 ± 1.0 

183.1 

181.0 

M°) 

4.4 ± 0.9 

2.3 

4.7 

v (km/s) 25.6 ± 0.8 

23.0 

28.9 

O 

356.6 

356.6 

356.6 


'Sekanina (1973), 2 Whipple (1957) 


Cluster of 5 seen on 17 Mar 2013 



Orbits of the cluster of 5 were very similar 
with the following average orbital elements: 



meteoroids 

NVI 

EVI 

a (AU) 

2.25 ±0.17 

1.69 

2.82 

e 

0.79 ± 0.02 

0.71 

0.86 

i(°) 

5.26 ± 1.02 

3.7 

5.2 

<o(°) 

280.32 ±2.11 

282.4 

285.8 

Q(°) 

356.65 ± 0.07 

358.0 

355.1 

q (AU) 

0.48 ± 0.02 

0.496 

0.40 

Q(AU) 

4.0 ±0.3 

2.89 

5.25 

Tj 

3.1 ±0.2 

Indicates -asteroidal body 
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Favorable Virginid 
radiant geometry 


^IflLlmpactFlash 


subradiant 






m 


ft 


Pink indicates the portion of the moon visible to the radiant. 
Impact angle -56° from horizontal. 
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Putting it all together 



All- sky meteor cameras detected a deeply 
penetrating cluster of 5 fireballs on 17 March. 

Radiant and orbital elements consistent with the 
Virginid Meteor Complex (EVI/NVI). 

Impact flash rate increased to 1 every 0.87 hours 
on 17 March. (4 impacts in 3.5 hours) 


Impact Constraints 

Assume impact flash was part of 
Virginid Meteor Complex 
v g = 25.6 km/s 
% = 56 ° 

^ Asteroidal? (Tj = 3.1 ± 0.2) 
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Mapping the impact location 

"Rough workflow" 



Clementine basemap 



ArcMap was used to georeference the 
lunar impact 3 times, at peak brightness 
and late impact. 


Using the geometric center of the flash 


Euler 

MARE IMBRIUM 

r. f 

' Late ^ 

.J impact 

- •*+* c 

Pytheas 


Rough workflow 
Avg late impact 

v. 

Draper 

50 km 

◄ ► 



Average predicted crater position 
20°.60 ± 0.17 N, 23°.92 ± 0.30 W 
was sent to LRO. 
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Mapping the impact location 

"Refined workflow" 



LRO basemap 



ArcMap was used to georeference the 
lunar impact following the geolocation 
workflow. 


Using the intensity-weighted center of the flash 


Euler 

Refrac corr. 

’ • \ 

MARE IMBRIUM 

V ‘ 

* ' 'S 

\ > ■ ) 

/ 

Pytheas 

Refined 


workflow 

0 12.5 ^5 r 50 Kilumu .uiy 

J 1 'i - ^ 1 i t.., ’ 

» * i ^ 

t « „ 

Draper 


Refined: Nominal predicted crater position 

20°. 6644 N, 24°. 1566 W 


Refrac corn Final predicted crater position 

20°. 6842±g;l||f N, 24°. 2277±g;l|g^ W 
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Impact crater found by LRO! 

Robinson et al. (2014) 


March 1 7 th Impact 



Image from Robinson (2013) 
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Features 

• Fresh, bright ejecta 

• Circular crater 

• Asymmetrical ray pattern 


Crater info 

• Rim-to-rim diameter = 18 m 

• Inner diameter = 15 m 

• Depth ~ 5 m 

Actual crater location 

• 20.7135°N, 24.3302°W 

Impact Constraints 
^ Circular crater, impact 
angle constrained 0 h >15° 
^ Ejecta gives no azimuth 

constraint (Robinson, personal comm.) 
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Comparison of geolocation 
results to obs crater location 






Method 

Longitude 
(° W) 

Latitude 

(°N) 

Angular distance 
from observed (°) 

Surface distance 
from observed (km) 

Rough workflow 

23.922 

20.599 

0.39875 

12.096 

Refined workflow 

24.1566 

20.6644 

0.169665 

5.1469 

Refined, with refraction correction 

24.2277 

20.6842 

0.100261 

3.0415 

LRO observed 

24.3302 

20.7135 

- 

- 
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Crater scaling laws 



Gault's scaling law (Gault 1974) for D < 100 m 
D = O.25p/ 167 p t _o - 5 E o - 29 (sin0 ft ) 1/3 

D = transient crater diameter 

p p = projectile density 

p t = target density 

E = kinetic energy of projectile 

6 h = impact angle measured wrt horizontal 


Holsapple crater calculator 

(http://keith.aa.washington.edu/craterdata/scaling/index.htm) 


Crater Sizes from Explosions or Impacts 

by 

About.. .=> I Keith A. Holsapple « 



Default values for above choice. Choose "Other" to change: 

Mass density: |i 5 1 1 gm.tm*3 ▼ | 

Target Cohesion (lab scale): |iex 1 1 a>nexm A 2 ▼ | 

Assumed Cohesion at this crater scale: 

Friction Angle: [33 1 1 degrees ▼ | 

Porosity (%): [to 1 1 percent ▼**[ 


The Impactor: 


I mpactor Material: 

I c-Type ▼ 1 



The Excavation Crater Size is: 


Crater Diameter 



Depth 

meters 

Rim Diameter 

^B^e^^Trneters^^^^l 

Lip Height 

^^g0 359 | meters 


Excavation Volume 


Ejecta Mass 

[5 03E- *9 rams T 


Ejecta Characteristics: 


The mass ejected with velocity greater than 
Is:bt percent of total, or E5* 
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Crater scaling laws 



• Transient crater diameter • Crater rim diameter 

• Measured at the height of the • Measured from rim to rim 

pre-existing surface • Aka outer diameter 

• Aka inner diameter, apparent 
diameter 

Holsapple assumes D rim = 1.3 D similar to Melosh (1989) 
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JAC P^ Transient crater diameter estimates 

Assumptions: Virginid v gfoc =25.7 km/s, 0 h = 56°; p t = 1500 kg/m 3 (regolith) 


Model 

Lum eff. r| 

KE 

xlO 9 (J) 

Mass 

(kg) 

Pp 

(kg/m 3 ) 

Dcaic 

(m) 

Dobs 

(m) 

% Err 

Gault's crater 
scaling law 

(Gault 1974) 

5 X 10' 4 

(Bouley et al. 2012) 

14 

[9.4,22] 

42 

[28,66] 

1800 

18.5 [16.5,21.1] 

15 

23% 

3000 

20.2 [18.0,23.0] 

15 

35% 

1.3X10- 3 

(Moser etal. 2011) 

5.4 

[3.6, 8.4] 

16 

[11,26] 

1800 

14.1 [12.5,16.0] 

15 

6% 

3000 

15.3 [13.6,17.4] 

15 

2% 

Holsapple's 

online 

calculator 

(Holsapple 1993) 

5 X 10" 4 

14 

[9.4,22] 

42 

[28,66] 

1800 

12.2 [10.9,13.8] 

15 

19% 

3000 

12.5 [11.1,14.2] 

15 

17% 

1.3X10 3 

5.4 

[3.6, 8.4] 

16 

[11,26] 

1800 

9.3 [8.3,10.5] 

15 

38% 

3000 

9.5 [8.5,10.8] 

15 

37% 


Two example values of r\ from the literature yield large ranges for KE and mass. 
Consequently, model results are highly dependent on luminous efficiency r|. 


Assuming a velocity dependent r\ = 1.3xl0 3 , these model 
results are consistent with the observed crater diameters. 

D ca i c = 8-18 m transient crater D obs = 15 m inner (‘transient’) 

D calc = 10-23 m rim-to-rim D obs = 18 m rim-to-rim 
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fySSJ Luminous efficiency: 



D t = 15 m from crater measurements "1 Gault’s crater scaling law (Gault 1974) rearranges to give 

LE A = 7. 1 x 10 6 J from flash measurements J jj a vs 0 h without assuming impact speed. 

nx = LE a /KE = LE A / ( 4.0 D t 0167 P f 0 - 5 sin - 1/3 ^) 1/0 - 29 



Projectile densities 

* pp = 1000 kg/m J 

— jt? p = 1 S00 kg/m J 

— pp = 2B00 kg/m 3 
pp = 3000 kg/m 3 

— pp= 3300 kg/m 3 


Typical values of rj A derived from lunar regolith range from 2 ± 1 x 1 O' 4 to 2 ± 1 x 10' 3 . 
Assuming association with the Virginids, 6 h = 56° and 7.52.5 x 10' 4 < iJ A < 1.52!* x 10' 3 . 
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Impact summary 


Date of impact: 
Duration of impact: 
Corrected flash peak R magnitude: 
Luminous energy generated by impact: 
Estimated kinetic energy of impactor: 
Estimated mass of impactor: 
Estimated diameter of impactor: 

Crater diameter: 
Crater location: 
Possible meteor shower association: 


17 March 2013 3:50:54 UTC 
1.03 s 

3.0 ±0.4 
7 . 1 S 4 X 10 6 J 

5.4-u x 10 9 J = 1.3 tons of TNT (assuming r| = 1.3 x 10" 3 ) 

16? kg (assuming v = 25.7 km/s) 

22 ± 3 cm (assuming p p = 3000 kg/m 3 ) 

18 m rim-to-rim, 15 m inner (‘transient’) 
20.7135° N, 24.3302° W 

Virginid Meteor Complex 
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Ejecta distribution 

after Robinson et al. (2014) 



Ejecta in multiple reflectance "zones" 

High reflectance zone 10-20 m SW, <10 m NE 
Low reflectance zone 50 m WSW, 80 m ENE 

-300 m rough semicircle 
Low reflectance zone -1 km centered in NE 
248 circ/irreg splotches within 30 km see Robinson etai ( 2014 ) 

& ± for more details 


Impact Constraints 

^ Circular crater, impact angle constrained >15° 


HRZ - impact possible from SE or NW 
LRZ - impact possible from SW 

.'. no azimuth constraint (Robinson, personal comm.) 


An impact from the SW is consistent with an 
impactor from the Virginid Meteor Complex. 
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